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SUMMARY 

Many  applied  problems  would  benefit  from  being  able  to  describe  postfis- 
sion  decay  of  beta  or  gamma  spectra  accurately  over  a limited  range  of  times 
by  only  a small  number  of  parameters.  This  report  describes  the  preparation 
of  such  a set  by  least-squares  fitting  of  a series  of  exponential  functions 


to  each  group  in  selected  multigroup  spectra  calculated  by  summation  of  fis- 
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sioo-product  data  in  ENDF/B.  The  parameters  are  for  fission  of  U,  U, 


and  Pu  induced  by  fission-spectrum  and  14-MeV  neutrons.  They  cover  the  time 
range  from  0.1  to  1000  seconds  after  fission. 

Approximately  45 % of  the  photons  emitted  from  zero  to  infinite  decay  time 
emerge  in  the  first  millisecond  after  fission.  These  are  not  included  in  the 
data  set  described  in  this  report.  Of  all  beta  particles  emitted  from  the 
instant  of  fission,  or  all  gamma  rays  emitted  after  the  first  millisecond, 
about  70%  emerge  during  the  0.1-to-1000-second  interval  covered  by  our  fits. 
Extrapolation  of  our  fits  to  infinite  time  misses  an  average  of  about  5%  of 
the  postmillisecond  yield. 

The  ENDF/B  data  have  uncertainties  in  total  energy-release  rate3  of  less 

than  10%  for  decay  times  greater  than  about  100  seconds.  The  fits  can  achieve 

few-percent  accuracy  with  five  or  fewer  terms  in  the  series.  At  the  expense 

of  one  additional  term,  a given  fitting  accuracy  can  be  achieved  with  a more 

compact  parameter  set  that  uses  a single  set  of  decay  constants  per  group  for 

all  six  combinations  of  fissioning  isotope  and  neutron  energy  but  individual 

amplitudes  for  each  combination.  We  find  the  earlier  work  by  Dieckhoner  on 

235 

the  beta  spectrum  from  U irradiated  by  fission-spectrum  neutrons  to  oe  in 
moderately  good  agreement  with  modern  data  near  the  peak  in  the  spectrum  for 
several  decay  times,  but  in  poor  agreement  elsewhere.  We  anticipate  a sub- 
stantial improvement  in  reliability  of  calculations  like  ours  when  Version  5 
of  ENDF/B  is  released,  particularly  for  the  short  decay  times  emphasized  in 
this  work. 
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PREFACE 

The  work  described  in  this  report  would  have  been  impossible  without  the 
preexisting  code  system  and  data  base  developed  with  the  support  of  several 
agencies  of  the  U.S.  Government.  The  fission-product  data  in  ENDF/3  are  the 
fruit  of  a long-term  program  of  the  Division  of  Reactor  Development  and  Dem- 
onstration of  the  U.S.  Energy  Researcn  and  Development  Agency  (ERDA)  (now  De- 
partment of  Energy).  The  U.S.  Nuclear  Regulatory  Commission  supported  the 
evolution  of  the  computer  program  CINDER-7  into  CINDER-10,  the  development  of 
the  programs  FPDCYS  and  FPSPEC,  and  the  use  of  CINDER- 10  in  processing  five 
of  the  six  basic  data  sets  used  in  the  present  work.  AFWL,  to  whom  this  report 

is  addressed,  supported  the  reduction  of  preliminary  Version-5  ENDF/B  data  to 
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individual-nuclide  fission  yields  for  14-MeV  neutrons  on  Pu  and  the  corre- 
sponding processing  through  CINDER.  AFWL  also  supported  processing  of  all  of 
the  data  through  a streamlined  version  of  FPSPEC,  as  well  as  the  development 
and  use  of  the  program  FPSPFT  and  its  ancillary  codes  (which  constituted  the 
main  part  of  the  present  work).  ERDA's  Division  of  Military  Applications 
furnished  supporting  services  in  the  preparation  of  this  report. 
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SECTION  I 
INTRODUCTION 


This  report  describes  the  preparation,  for  the  Air  Force  Weapons  Labora- 
tory (AFWL) , of  a compact  parameterization  of  the  beta  and  gamma  spectra  emit- 
ted by  fission  products  during  the  interval  from  0.1  to  1000  seconds  after  a 
very  short  fission  burst.  Although  this  project  was  designed  specifically  to 
supply  input  to  computer  codes  that  calculate  the  injection  of  electrons  from 
a high-altitude  nuclear  explosion  into  the  earth's  magnetosphere,  the  results 
have  obvious  applications  to  other  problems.  Indeed,  the  gamma  spectra  were 
included  without  any  immediate  objective  in  mind,  taking  advantage  of  the  fact 
that  they  constituted  an  inexpensive  add-on  to  the  beta  spectra.  The  param- 
eterization that  we  have  used  is  least-squares  fitting  of  a series  of  expo- 
nential functions. 

The  specific  task  for  which  this  work  was  performed  has  had  the  effect  of 

limiting  its  scope  in  two  respects.  In  the  first  place,  since  it  is  intended 

to  be  applied  to  nuclear  explosions,  the  resulting  parameter  set  is  restricted 
235  238  239 

to  fission  of  U,  U,  and  Pu  by  fission  neutrons  and  by  14-MeV  neutrons. 
Although  data  are  available  for  thermal  neutrons,  these  have  not  been  included 
In  the  set  because  practical  nuclear  explosives  make  negligible  use  of  thermal 
neutrons.  In  the  second  place,  the  time  range  of  0.1  to  1000  seconds  after 
fission  excludes  interesting  phenomena  at  both  earlier  and  later  times.  Both 
limits  were  chosen  to  permit  a quick,  accurate  solution  with  a minimum  number 
of  parameters.  Thus,  although  the  electron  spectrum  should  extrapolate  accu- 
rately to  zero  time,  the  parameterized  photon  spectrum  does  not  include  the 
intense  prompt  burst  or  succeeding  decay  of  metastable  states  of  the  primary 
fission  products  prior  to  one  millisecond,  which  constitute  the  source  of  the 

electromagnetic  pulre  that  follows  a nuclear  explosion.  Both  electron  and 

3 4 

photon  yield  curves  exhibit  a pronounced  kink  between  10  and  10  seconds 
after  fission.  Fitting  these  would  have  required  several  additional  parame- 
ters in  order  to  describe  fewer  than  15%  of  the  electrons  injected  into  the 
magnetosphere,  and  would  have  increased  the  cost  and  duration  of  this  project 
in  exchange  for  minimal  practical  gain. 

Our  basic  approach  is  usually  described  as  a summation  calculation. 

It  continues  the  trend  exemplified  by  the  earlier  work  of  Stovall  (ref.  1) 
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and  of  Dieckhoner  (ref.  2).  Rather  than  work  from  direct  measurements  of  time- 
dependent  spectra,  of  which  there  are  not  nearly  enough  to  cover  the  require- 
ments of  this  project,  we  have  reconstructed  the  spectra  from  detailed  data 
on  each  step  in  the  fission-and-decay  process.  This  was  possible  within  the 
existing  budgetary  constraints  only  because  other  sponsors*  had  already  sup- 
ported most  of  the  data  preparation,  and  the  prepared  data  were  stored  in  the 
Central  Computing  Facility  of  the  Los  Alamos  Scientific  Laboratory  (LASL)  in 

machine-readable  form.  The  only  basic  data  added  specifically  for  this  pro- 

239 

ject  were  fission  yields  for  14-MeV  neutrons  irradiating  Pu. 

Section  II  of  this  report  describes  this  data  base  in  greater  detail,  and 
summarizes  the  initial  computational  procedure,  including  the  parts  that  had 
already  been  completed  before  this  project  began.  Section  III  gives  a general 
overview  of  the  resulting  time-dependent  spectra.  Section  IV  discusses  the 
extensive  program  of  validating  the  data  and  the  computational  procedures  by 
comparison  to  benchmark  experiments.  Section  V gives  an  explicit  comparison 
of  the  present  results  with  those  of  Dieckhoner  (ref.  2).  Section  VI  de- 
scribes the  form  and  procedure  for  parameterizing  the  time-dependent  spectra 
for  use  by  AFWL.  Section  VII  discusses  extensions  and  improvements  to  both 
the  data  and  the  fitting  process  that  could  profitably  be  exploited  in  the 
future.  As  an  aid  to  actual  use  of  the  AFWL  parameter  set,  we  have  included 
appendices  that  describe  the  mathematics  of  the  fit  and  the  format  of  both  the 
output  listings  from  the  fitting  program  and  the  cards  used  for  transmitting 
the  fitted  parameters. 

Since  we  shall  refer  frequently  to  energy-isotope  cases;  that  is,  to 
specific  combinations  of  neutron  energy  and  fissioning  nuclide,  we  shall  hence- 
forth use  the  abbreviations  for  these  cases  listed  in  table  1.  The  last  column 
of  the  fable  lists  the  number  of  individual  decay  times  for  which  we  have  spec- 
tra available  for  fitting  in  the  time  range  between  0.1  and  1000  seconds  after 
fission. 


*See  the  preface  to  this  report  for  the  contributions  supported  by  each 
sponsor. 
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f Table  1 

| abbreviations  used  to  designate  data  sets 


Fissioning 

Nuclide 

Neutron 

Energy 

Abbreviation 

Number  of 
Decay  Times 
0.1-1000  8 

235 

thermal 

25T 

29 

239_ 

Pu 

thermal 

49T 

29 

235u 

fiss*  spec* 

25F 

9 

238u 

fiss.  spec.* 

28F 

29 

239Pu 

flss.  spec.* 

49F 

29 

235u 

14  MeV 

25H 

9 

238u 

14  MeV 

28H 

9 

239Pu 

14  MeV 

49H 

29 

*Fl88.  spec.  • fission-spectrum  neutrons 
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SECTION  II 

CALCULATION  OF  TIME-DEPENDENT  SPECTRA 

All  of  the  data  used  in  this  work  are  taken  from  the  Evaluated  Nuclear 

Data  File,  part  B (ref.  3),  which  is  customarily  abbreviated  ENDF/B.  The 

239 

individual-nuclide  yields  of  the  fission  products  from  14-MeV  fission  of  Pu 
(case  49H)  are  taken  from  a preliminary  set  for  Version  5 of  ENDF/B;  all  other 
data  are  from  Version  4.  Table  2 summarizes  the  nature  and  extent  of  the  data 
in  ENDF/B-4.  Most  of  the  data  are  from  actual  measurements,  but  some  of  the 
half-lives  and  decay  energies  have  been  taken  from  nuclear  systematics,  espe- 
cially for  nuclides  with  half-lives  less  than  one  second.  In  Version  5 im- 
proved models  have  bean  used  to  supplement  measurements  in  distributing  the 
mass-chain  yields  among  the  individual  nuclides  In  each  chain.  In  particular, 
Madland  and  England  (ref.  4),  find  that  a correction  for  the  effect  of  the 
nuclear  pairing  force  varies  greatly  from  one  fissioning  nuclide  to  another 
for  low  neutron  energies,  but  has  only  a minor  effect  at  14  MeV.  In  ref.  5 
the  same  authors  describe  a semiempirical  model  for  the  ratio  of  yields  to 
the  ground  and  isomeric  states  of  the  same  nucleus.  In  the  past  it  has  been 
assumed,  in  the  absence  of  measurements,  that  the  two  states  are  populated 
equally,  whereas  Madland  and  England  find  that  the  isomeric  state  is  usually 
strongly  favored. 

Table  2 

SUMMARY  OF  ENDF/B-4  FISSION-PRODUCT  DATA 

825  nuclides  (total,  counting  isomers  separately) 

112  isomers  with  half  lives  > 0.1  second 
42  different  elements 
96  different  mass  numbers 

77  nuclides  marked  as  gaseous  isotopes  of  4 elements 

181  have  neutron- interaction  cross  sections 

180  have  explicit  data  on  beta  and/or  gamma  spectra 

712  are  unstable  and  have  average  a,  $,  and  y energy 
and  branching  fraction 

825  have  fission  yields  for  each  of  6 fissionable 
nuclides  (~  10^  yields)  for  one  or  more  ranges 
of  neutron  energy 

10 


AFWL-TR-78-4 


ENDF/B  data  were  originally  strongly  oriented  towards  the  needs  of  the 
reactor-design  community.  Accordingly,  experimental  data  and  theoretical 
interpolations  now  in  the  files  emphasize  energy  production  rather  than  par- 
ticle spectra.  The  spectral  data  are  concentrated  on  times  greater  than  a 
minute  after  fission  and  on  fission  induced  by  thermal  neutrons.  Thus,  the 
data  do  not  serve  our  particular  needs  as  well  as  they  do  the  needs  of  reactor 
design. 

The  primary  code  that  LASL  uses  (ref.  6)  for  generating  the  time  depend- 
ence of  fission-product  decay  is  CINDER  (currently  CINDER- 10).  Although  it 
was  developed  earlier  and  entirely  independently,  CINDER  contains  substantially 
all  of  the  features  recommended  by  Dieckhoner  (ref.  2),  in  addition  to  having 
access  to  the  vast  amount  of  recent  data  in  ENDF/B.  It  calculates  production 
and  depletion  of  fission  products  by  applying  a known  neutron-flux  history  to 
multigroup  cross  sections  taken  from  ENDF/B.  In  the  present  application  to 
essentially  instantaneous  fission  we  have  uniformly  used  a constant  flux  of 
10J"3/(s*cm2)  applied  for  10  seconds.  Decay  is  followed  along  the  mass 
chains,  as  well  as  between  chains  when  delayed  neutrons  are  emitted.  The 
actual  coupled  decay  scheme  of  825  nuclides  is  decomposed  into  linearized 
chains  containing  2811  nuclides,  of  which  1722  are  out-of-sum  virtual  mem- 
bers to  eliminate  double  counting  of  partial  yields  and  the  rest  are  summed 
to  give  the  actual  yields.  CINDER  solves  the  linearized  chains  using  a com- 
pletely general  solution  to  the  coupled  differential  equations.  Because  some 
of  the  chains  contain  as  many  as  twenty  members,  very  careful  attention  has 
been  given  to  preventing  the  accumulation  of  rounding  errors.  CINDER  output 
lists  the  concentration,  activity,  and  energy-emission  rate  for  each  nuclide 
(Including  isomers)  separately  as  well  as  their  overall  sums.  In  addition, 
gaseous  and  volatile-solid  isotopes  are  summed  separately. 

The  180  nuclides  that  have  detailed  spectra  in  ENDF/B  have  been  processed 
separately  by  the  code  FPDCY3  (ref.  7).  FPDCYS  uses  the  exact  Fermi  express- 
ions to  generate  beta  spectra  from  the  tabulated  end-point  energies,  transi- 
tion types,  and  branching  ratios,  then  calculates  the  integrals  over  75  100-keV 
bins.  The  tabulated  gamma-ray  line  spectra  from  ENDF/B  are  simply  binned  into 
125  50-keV  groups.  Both  sets  of  spectra  are  stored  for  use  by  other  codes. 

For  our  project  we  have  sorted  the  spectra  into  a standard  order  and  placed 
them  in  a compact  binary  library  for  rapid  retrieval.  Similarly,  we  have 
picked  up  from  the  stored  output  of  CINDER  the  individual-nuclide  data  as  a 
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function  of  decay  time,  sorted  them  into  the  same  standard  sequence,  and  stored 
them  in  another  binary  library  along  with  the  parameters  that  describe  the 
irradiation  history. 

The  final  step  in  generating  the  composite  spectra  is  to  combine  the 
stored  spectra  from  FPDCYS  with  the  stored  irradiation  parameters  and  time- 
dependent  decay-rates  from  CINDER,  using  a code  named  FPSPEC  (ref.  7).  The 
result  is  absolute  particle-emission  rates  per  fission  from  the  180  nuclides 
with  known  spectra  (ref.  8).  The  CINDER  calculations  also  supply  the  total 
electron  and  photon  energy-emission  rates  for  these  nuclides,  and  independently 
for  the  entire  ensemble  of  825  fission  products.  Since  we  have  no  further 
information  available,  we  assume  that  the  average  spectrum  of  the  remaining 
645  nuclides  is  the  same  as  that  of  the  180,  and  normalise  the  calculated 
spectrum  to  the  total  beta  or  gamma  energy-emission  rate.  We  shall  consider 

the  validity  of  this  normalization  in  Sections  III  and  IV.  FPSPEC  stores 
these  normalized  spectra  in  a cumulative  binary  library  of  its  own,  to  serve 
as  input  to  the  exponential-fitting  process. 

The  CINDER  output  used  in  this  project  was  generated  at  various  times 
during  1S76  and  1977  for  various  immediate  purposes.  The  earliest  data  sets 
contained  calculations  only  at  decay  times  of  1 and  5 times  powers  of  10 
seconds.  Using  the  notation  of  table  1,  sets  25F,  25H,  and  28H  are  of  this 

3 

type.  These  "sparse"  sets  have  only  nine  times  in  the  range  0.1  to  10  seconds, 
which  seriously  restricts  the  fitting  operation.  More  recently,  sets  28H  and 
49F  were  rerun  to  add  points  at  1.5,  2,  4,  6,  and  8 in  each  time  decade,  rais- 
ing the  total  to  29  times  per  set.  We  shall  call  these  "dense"  sets.  The  only 
case  run  expressly  for  this  project,  49H,  is  also  dense,  but  in  each  decade  we 
have  substituted  a point  at  3 for  the  one  at  5 to  produce  a more  even  loga- 
rithmic spacing.  The  number  of  decay  times  per  case  is  included  in  table  1. 
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be  calculated  by  CINDER.  The  shape  of  the  spectra,  on  the  other  hand,  is 
reasonably  certain  only  in  proportion  to  the  fraction  of  the  total  energy  re- 
lease that  occurs  in  the  180  nuclides  for  which  ENDF/B-4  contains  detailed 
spectral  data.  Figures  1 and  2 show  this  fraction,  for  electrons  and  photons, 
respectively,  for  each  of  our  six  cases.  The  lower  part  of  each  figure  ap- 
plies to  fission  induced  by  a fission-neutron  spectrum  (cases  25F,  28F,  and 
49F) , The  upper  part  displays  the  ratio  of  the  14-MeV  fraction  to  the  cor- 
responding fission-energy  fraction  (25H/25F,  28H/28F,  or  49H/49F).  We  see 
that  in  general  less  than  20%  of  the  total  energy  is  in  known  spectra  at  0.1 
second  after  fission  (the  lowest  fraction  is  less  than  3%,  for  photons  in  case 
28F).  The  known  fraction  rises  slowly  until  about  20  seconds  after  fission 
for  photons  or  - 40  seconds  for  electrons,  then  rises  more  rapidly  until  about 
100C  seconds,  after  which  it  levels  off  well  above  90%  for  all  cases. 

It  is  tempting  to  assume  that  all  properties  of  fission-product  decay  at 
early  times  will  vary  slowly  and  smoothly  with  the  mass  of  the  fissioning 
nuclide  and  the  energy  of  the  incident  neutrons,  simply  because  the  number 
of  fission  products  is  very  large.  Figures  1 and  2 are  our  first  indication 
that  this  is  not  true,  and  we  shall  see  in  the  rest  of  this  report  that  var- 
iations of  up  to  a factor  of  two  between  cases  are  typical  of  most  quantities 
of  interest. 

Even  so,  we  can  make  two  generalizations  about  the  known-spectrum  frac- 
tion. First,  particularly  for  beta  decay  (fig.  1),  this  fraction  is  usually 

smaller  at  all  decay  times  for  14-MeV  neutrons  than  for  fission-spectrum  neu- 

238 

trons.  Secondly,  at  early  times  the  data  for  U are  conspicuously  deficient * 
2Y 

whereas  U is  the  first  to  become  well-known.  Both  tendencies  are  obvious 
consequences  of  the  theraa.1! -reactor  orientation  of  the  ENDF  system. 

Figures  3 and  4 show  the  total  particle  multiplicity  as  a function  of 
decay  time  for  electrons  and  photons,  respectively.  One  extreme  case  from 
each  figure  is  reproduced  in  the  other  to  emphasize  differences  between  beta 
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Figure  1.  Fraction  of  Total  Electron  Energy,  as  a Funci 
Time,  that  is  Emitted  by  the  180  Nuclides  wii 
Data  points  are  not  shown  after  1000  seconds, 
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Fraction  of  Energy  in  Known  Spectra 
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Figure  2.  Fraction  of  Total  Photon  Energy,  as  a Function  of  Decay  Time 
that  is  Emitted  by  the  180  Nuclides  with  Known  Spectra.  Dat 
points  are  not  shown  after  1000  seconds. 
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Figure  3.  Total  Intensity  of  Electrons  as  a Function  of  Decay  Time 
One  curve  for  photon  decay  is  included  for  comparison. 
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one  photon  per  electron  in  all  cases,  but  the  relative  difference  decreases 

steadily  for  the  first  hundred  seconds.  The  total  yield  of  either  particle 

is  systematically  of  the  order  of  10%  less  at  14  MeV  than  for  fission  neutrons. 

238  239 

In  both  figures  the  early  yield  from  U is  the  largest  and  from  Pu  the 
smallest,  but  the  difference  is  small  after  100  seconds.  The  decay  rate  for 
case  49F  is  markedly  smaller  than  for  the  other  five  cases. 

In  section  II  we  pointed  out  that  the  output  library  from  FPSPEC  contains 
absolute  emission  spectra  in  100-keV  bins  for  electrons  and  50-keV  bins  for 
photons.  In  order  to  reduce  this  information  to  manageable  size,  we  have 
rebinned  these  spectra  into  group  structures  appropriate  to  AFWL's  specific 
needs.  Table  3 lists  the  three  structures  adopted  for  this  purpose.  They 

Table  3 

AFWL  GROUP  STRUCTURES  FOR  FISSION-PRODUCT  SPECTRA 


Energy, 

38-gp 

17-gp 

20-gp 

Energy , 

38-gp 

17-gp 

20-gp 

MeV 

gamma 

gamma 

beta 

MeV 

gamma 

gamma 

beta 

7.5 

1 

1.00 

19 

11 

14 

7.0 

2 

0.95 

20 

6.5 

3 

.90 

21 

6.25 

i 

i 

.85 

22 

6.0 

4 

.80 

23 

5.5 

2 

2 

5 

.75 

24 

12 

15 

5.0 

3 

3 

6 

0.70 

25 

4.5 

4 

4 

7 

.65 

26 

4.0 

5 

5 

8 

.60 

27 

3.5 

6 

6 

9 

.55 

28 

3.0 

7 

7 

10 

.50 

79 

13 

16 

2.5 

8 

8 

11 

0.45 

30 

2.0 

9 

9 

12 

.40 

31 

1.9 

10 

.35 

32 

1.8 

11 

.30 

33 

1.7 

12 

.25 

34 

14 

17 
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14 

10 

13 

0.20 

35 

1.4 

15 

.15 

36 
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16 

■ 

.10 
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consist  of  a coarse-group  structure  for  both  electrons  and  photons  and  a 
compatible  second  fine-group  structure  for  photons.  We  shall  use  these  group 
structures  to  display  a sampling  of  the  behavior  of  the  spectra. 

Figure  5 shows  20-group  electron  spectra  for  the  25F  case,  taken  at  even 

4 

powers  to  ten  from  0.1  to  10  seconds  after  fission.  The  early  spectrum  has  a 
maximum  near  2 MeV  that  is  roughly  a factor  of  5 greater  than  the  intensity  in 
the  lowest-energy  bin.  The  spectrum  cuts  off  arbitraily  at  7.5  MeV,  where  the 
original  binning  for  the  180  known  spectra  ends  in  FPDCYS  (ref.  7).  The  inten- 
sity at  the  cutoff  energy  is  approximately  1%  of  the  maximum  at  0.1  second  and 
decreases  rapidly  at  later  times.  At  later  times  the  maximum  moves  towards 
lower  energies  and  the  high-energy  part  of  the  spectrum  drops  more  steeply. 
After  10  seconds  the  peak  disappears  altogether,  leaving  only  the  monotonic 
decrease  with  energy. 

The  corresponding  38-group  photon  spectra  for  25F  appear  in  fig.  6.  Here 
the  individual  lines  produce  such  large  fluctuations  that  it  is  difficult  to 
plot  them  in  a single  figure  even  at  time  intervals  of  a factor  of  ten.  In- 
deed, we  have  had  to  omit  the  spectrum  at  one  second  after  fission.  The  aver- 
age energy  of  the  spectrum  at  early  times  is  conspicuously  lower  than  for 
electron  emission,  with  the  four  highest  peaks  all  below  1 MeV  and  a sharp 
drop  in  intensity  above  1.5  MeV.  The  early  intensity  in  the  lowest-energy 
bin  is  a factor  of  a thousand  lower  than  in  the  first  peak,  but  this  extreme 
difference  disappears  after  ten  seconds.  By  100  seconds  after  fission  a broad 
background  is  evident  under  much-reduced  fine  structure,  with  a maximum  near 
0.6  MeV.  The  high-energy  tail  of  the  spectrum  is  quite  flat  at  early  times, 
but,  as  in  the  electron  spectrum,  becomes  much  steeper  at  later  times.  Note 
that  ENDF/B  does  not  list  any  photons  more  energetic  than  6.25  MeV. 

The  densities  of  many  individual  fission-product  nuclides  exhibit  growth 
followed  by  decay,  but  in  the  electron  spectrum  summed  over  all  the  nuclides 
none  of  the  100-keV  bins  exhibits  growth.  On  the  other  hand,  because  they 
are  easily  dominated  by  one  nuclide  many  of  the  50-keV  photon  bine  do  exhibit 
a maximum  in  the  decay  curve.  Figure  7 shows  an  example;  namely,  the  30-to- 
100  keV  bln  for  case  49F.  The  same  bin  in  the  electron  spectrum  has  no  cor- 
responding maximum.  In  contrast  to  the  total  intensity  as  shown  in  figs.  3 
and  4,  in  this  bln  the  photon  intensity  actually  drops  below  the  electron 
Intensity  after  600  seconds. 
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Figure  5.  Selected  20-Group  Electron  Spectre  Following  Fission  of  U 
by  Fission  Neutrons 
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Figure  6.  Selected  38-Group  Photon  Spectra  Following  Fission  of  235U 

by  Fission  Neutrons.  The  "extrapolated"  values  are  discussed 
in  the  text. 
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Figure  7.  Example  of  Growth  Followed  by  Decay  for  Low-Energy  Photons 
Following  Fission  of  ^39pu  by  Fission  Neutrons.  Note  that 
growth  is  not  observed  for  electrons  of  the  same  energy. 
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As  a final  example  of  diverse  behavior,  let  us  examine  the  0.5-to-0.55 
MeV  bin  in  the  photon  spectrum,  which  is  in  the  center  of  a peak  in  the  0.1- 
second  spectrum,  and  compare  it  to  the  broader  0.5-to-0.75  MeV  bin  in  the 
electron  spectrum.  In  figs.  8 and  9 we  see  that  the  photon  multiplicity  is 
roughly  thirty-fold  greater  than  the  electron  multiplicity  at  early  times,  but 
decreases  very  rapidly  and  crosses  the  electron  decay  curve  twice  in  the  first 
thousand  seconds.  The  slower  decay  in  case  &9F, which  we  pointed  out  in  figs. 

3 and  4,  is  especially  marked  in  fig.  9.  These  two  figures  show  only  the 
fission-energy  cases.  The  ratios  of  the  intensities  for  14-MeV  neutrons  to 
those  for  fission  neutrons  are  shown  in  the  two  parts  of  fig.  10.  The  vari- 
ations in  the  electron  intensities  are  comparable  to  those  that  we  observe  in 
the  total  intensity.  The  14-MeV/ fission  ratio  for  the  narrow  photon  bin,  on 
the  other  hand,  shows  large  fluctuations  for  all  three  fissioning  nuclei 
(of  the  order  of  30%) . 
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Figure  8.  Decay  Curves  for  a Single  Electron  Group  Following  Fission 
Induced  by  Fission  Neutrons.  The  dashed  line  is  che  corre- 
sponding curve  for  photons  from  ^39pu  in  the  center  of  the 
peak  at  0.5  Mev. 
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Figure  10*  Ratio*  of  the  Decay  Curve*  for  Fission  Induced  by  14-HeV 
Neutrons  to  those  Induced  by  Fission  Neutrons  Shoivn  In 
Figs.  8 and  9 
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SECTION  IV 


VALIDATION  OF  DATA 


As  with  the  data  themselves,  validation  of  ENDF/B-4  fission-product  files 

235 

has  emphasized  heat  production  following  U fission  induced  by  thermal  neu- 
trons (ref.  9).  A number  of  spectrum  benchmarks  exist  (discussed  in  ref.  10), 
but  all  of  them  are  likewise  for  thermal  neutrons.  Nevertheless,  it  is  pos- 
sible to  make  some  estimates  for  other  fissioning  nuclides  and  other  energies. 

Schmittroth  and  Schenter  (ref,  11)  have  propagated  the  experimental  and 

modeling  uncertainties  underlying  the  ENDF/B  data  through  the  summation 

calculations  in  order  to  estimate  the  a 'priori  uncertainties  in  the  total 

energy-emission  rate.  Table  4 shows  part  of  their  results  (excerpted  from 

235 

ref.  12)  for  thermal  neutrons  incident  on  U (case  25T) , as  well  as  for  our 
cases  28F  and  49F.  The  same  authors  (ref.  13)  have  performed  a least-squares 
analysis  of  the  decay-heat  benchmark  measurements  in  order  to  derive  a poster- 
iori estimates  of  the  same  uncertainties,  and  we  have  included  their  results 
for  25T  in  table  4.  The  dominant  uncertainties  in  the  ENDF/B  data  come  from 
the  nuclide-yield  model,  and  at  early  times  from  the  mass-excess  estimates 
for  nuclides  far  from  stability.  We  see  from  the  table  that  the  a priori 
uncertainties  are  similar  for  all  cases.  However,  for  25T  the  a posteriori 
uncertainties  are  dramatically  smaller  for  decay  times  between  2 and  106  sec- 
onds. Indeed,  they  turn  out  to  be  dominated  by  systematic  disagreements  among 
the  thermal  benchmarks. 

The  thermal-neutron  benchmarks  include  a beta-spectrum  measurement  made 
at  the  University  of  Illinois  (ref.  14)  following  a 15-ms  puJ.se  and  a gamma- 
spectrum  measurement  made  at  the  Oak  Ridge  National  Laboratory  (ref.  15)  fol- 
lowing a 1-s  pulse,  both  of  which  constitute  stringent  tests  of  the  short-half- 
life data.  The  Illinois  measurement  agrees  very  well  with  ENDF/B-4  after  100 
seconds.  At  shorter  decay  times  and  electron  energies  greater  than  0.5  MeV, 
beta-spectrum  measurements  from  the  two  laboratories  disagree  (ref.  12)  with 
each  other  by  amounts  comparable  to  their  disagreement  with  the  summation  cal- 
culations. The  Oak  Ridge  gamma-spectrum  measurements  are  in  good  agreement 
with  our  calculations  after  50  seconds,  except  for  a narrow  peak  near  0.4  MeV, 
which  we  systematically  underestimate  by  as  much  as  a factor  of  two  between 
20  and  200  seconds  after  fission. 
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Table  4 

PERCEN-  UNCERTAINTIES  IN  ENERGY-EMISSION  RATE  FOLLOWING  A FISSION  BURST* 


Decay 
Time,  s 

235Uf 

Thermal  neutrons 

238u, 

Fission  neutrons 

Fission  neutrons 

Est.** 

I Meas.*** 

Est.** 

Est.** 

0.0 

24.0 

32.1 

42.5 

0.1 

23.2 

31.4 

42.3 

0.2 

21.9 

28.2 

38.2 

0.5 

20.1 

24.7 

32.2 

1.0 

18.6 

28. 

22.4 

27.2 

2.0 

17.4 

8. 

20.7 

22.5 

5.0 

16.2 

19.1 

18.8 

10.0 

15.1 

2.5 

17.7 

17.1 

20.0 

13.6 

2.1 

16.3 

15.6 

50.0 

10.1 

13.4 

12.8 

1 x 102 

6.8 

1.7 

9.6 

9.8 

2 x 102 

5.1 

6.6 

7.2 

5 x 102 

4.8 

5.3 

5.8 

1 x 103 

4.6 

1.7 

4.7 

5.5 

2 x 103 

4.4 

4.5 

5.6 

5 x 103 

4.3 

4.5 

6.2 

1 x 10  4 

4.3 

1.6 

4.7 

5.9 

2 x 104 

3.7 

4.1 

4.6 

5 x 104 

3.2 

3.5 

3.9 

1 x 105 

3.2 

1.7 

3.3 

3.7 

2 x 105 

2.3 

2.4 

3.0 

5 x 105 

1.7 

2.0 

2.4 

1 x 106 

1.5 

2.0 

1.8 

2.0 

1 x 107 

1.6 

2.0 

2.1 

1.9 

1 x 108 

2.2 

2.0 

3.8 

3.8 

* Excerpted  from 

ref.  12. 

**  Uncertainty  estimated  by  propagating  ENDF/B  uncertainties  to  derived 
emission  rate. 


***  Uncertainty  deduced  from  least-squares  analysis  of  benchmark  measure- 
ments. 
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From  the  above  considerations,  we  estimate  that  our  calculated  spectra 
are  about  as  accurate  as  present-day  direct  measurements  for  decay  times  for 
which  more  than  half  of  the  yield  comes  from  the  180  nuclides  for  which  we 
have  spectra  in  ENDF/B.  From  figs.  1 and  2 we  see  that  this  threshold  varies 
between  50  and  100  seconds  after  fission.  The  accuracy  of  our  detailed  spectra 
at  earlier  times  depends  on  the  postulated  similarity  of  the  overall  spectrum 
to  that  of  the  180  nuclides,  so  that  for  the  present  we  know  only  that  it  is 
no  better  than  the  examples  shown  in  reference  10  (see  particularly  their 
figs.  52  to  62  and  91  to  94).  Our  total  energy-emission  rates  are  probably 
correct  to  substantially  better  than  5%  for  decay  times  greater  than  10  sec- 
onds. 
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SECTION  V 

COMPARISON  WITH  EARLIER  CALCULATIONS 

For  many  years  the  calculations  begun  by  Stovall  (ref.  1)  and  improved 
by  Dieckhoner  (ref.  2)  have  been  the  standard  reference  for  fission-product 
spectra  cited  in  calculations  for  the  Air  Force.  We  have  pointed  out  in  Sec- 
tion I that  the  LASL  calculations  are  basically  a modern  incarnation  of  the 
same  computational  approach.  The  LASL  fission-product  work  includes  a number 
of  improvements  in  addition  to  incorporating  the  results  of  post-1963  experi- 
mental work. 

Dieckhoner  used  83  mass-chain  yields  taken  from  smooth  curves,  whereas 
ENDF/B  uses  actual  measurements  for  95  mass  chains  with  all  Irregularities 
included.  ENDF/B  uses  an  empirical  model  similar  to  the  equal-charge-dis- 
placement hypothesis  that  he  used  for  calculating  unmeasured  direct  yields. 

In  Version  5 (which  we  have  used  only  for  case  49H),  ENDF/B  adopts  two  new 
models  (refs.  4 and  5),  one  to  include  pairing  effects  for  individual  yields 
within  a given  mass  chain,  and  another  to  describe  the  branching  ratio  between 
ground  and  isomeric  states  (which  normally  strongly  favors  the  isomer,  contrary 
to  the  previous  common  assumption).  Dieckhoner' s longest  linearized  chains 
contained  only  10  members,  whereas  CINDER  uses  up  to  20  and  includes  inter- 
chain transitions  that  result  from  neutron  capture  and  emission.  Also,  CINDER 
uses  a single  general  solution  to  the  linearized-chain  differential  equations, 
rarher  than  special  forms  for  each  length  of  chain.  We  have  also  given  ex- 
haustive attention  to  the  problem  of  numerical  stability  of  the  solutions, 
which  Dieckhoner  does  not  seem  to  have  addressed  at  all,  but  which  we  now 
know  can  cause  order-of-magr.itude  errors.  Furthermore,  instead  of  estimating 
the  spectrum  from  average  ratios  of  beta  end-point  energy  to  total  decay 
energy,  we  have  used  exact  spectra  of  180  dominant  nuclides  as  models.  Finally, 
we  have  added  the  data  required  to  generate  photon  decay  curves  to  those  re- 
quired for  electrons. 

CINDER  includes  the  detailed  treatment  of  irradiation  history  recommended 
by  Dieckhoner,  so  that  it  can  follow  neutron-induced  changes  in  fission-product 
concentrations.  However,  in  the  absence  of  such  transmutation  it  is  not  neces- 
sary to  have  this  ability  in  order  to  correct  for  either  finite  irradiation 
time  or  finite  counting  time.  LaBauv^  et  aZ.,  point  out  in  ref.  16  that  it 
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should  be  possible  to  use  the  results  of  a fit  such  as  we  shall  describe  in 
Section  VI  to  transform  the  results  of  a short-burst  calculation  analytically 
into  the  result  for  any  specified  finite  irradiation  and  counting  times. 

Only  one  of  Dieckhoner's  four  cases  corresponds  directly  to  one  of  ours; 
namely,  the  electron  spectrum  for  case  25F.  Two  of  his  other  cases  are  for 
thermal  neutrons,  which  have  negligible  effect  in  practical  nuclear  explosions, 
and  the  fourth  is  a "thermonuclear"  case  which  contains  an  unspecified  mixture 
of  isotopes  and  neutron  energies. 

Figure  11  compares  Dieckhoner’s  25F  case,  with  our  calculations,  which 
have  been  interpolated  to  his  decay  times  by  using  the  fits  described  in 
Section  VI.  As  one  would  expect  from  the  paucity  of  measurements  for  short 
half-lives,  his  calculations  have  missed  about  half  of  the  intensity  extrap- 
olated back  to  zero  time.  He  also  had  serious  difficulty  with  the  high-energy 
part  of  the  spectrum,  which  at  7 MeV  he  underestimated  by  an  order  of  magnitude 
for  0 through  20  seconds  and  overestimated  by  factors  of  10  to  40  between  60 
and  300  seconds.  On  the  other  hand,  his  calculations  at  2,  6,  and  10  seconds 
agree  with  ours  within  about  10%  from  zero  to  the  maximum  in  the  curve,  and 
the  2-second  curve  is  still  less  than  25%  low  at  5 MeV.  After  20  seconds  his 
calculations  are  about  a factor  of  2 low  at  zero  energy,  but  all  of  them  cross 
our  curves  near  4 MeV  and  end  up  much  higher  at  7 MeV.  The  300-second  curve 
agrees  with  ours  within  about  15%  from  0.5  to  4 MeV.  At  1 MeV,  all  of  his 
curves  except  the  one  at  0 seconds  agree  with  ours  within  30%,  but  5 out  of 
these  7 are  low. 
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SECTION  VI 

PREPARATION  OF  COMPACT  PARAMETER  SET 

1.  FITTING  ALGORITHM 

One  obvious  method  for  representing  a multigroup  spectrum  of  fission- 
product  radiation  as  a function  of  time  is  to  perform  a least-squares  fit  of 
each  group  to  an  exponential  series  of  the  form 

m .t 

it  _ V 'r.  _ 8l 


where  y is  the  intensity  in  group  g at  decay-time  t and  m is  the  number  of 

O 

terms  in  the  series.  Previous  work  at  LASL  (ref.  16)  has  shown  that  this 
gives  a satisfactory  description  for  the  few-group  energy-emission  rate  and 
the  total  energy-emission  rate.  Indeed,  a 23-term  fit  to  the  total  energy 
has  been  proposed  (refs.  12  and  17)  as  an  engineering  standard  for  reactor 
use.  Our  purpose  in  the  present  project  was  to  extend  this  technique  to  the 
particle  spectra,  using  finer  group  structures  specified  by  AFWL. 

Our  basic  tool  for  the  fitting  was  a Fortran  program  called  FPSPFT  that 
was  written  specifically  for  this  purpose  and  run  on  a CDC-7600  computer  under 
the  Chili  Ridge  Operating  System  (CROS) . The  input  data  were  prepared  under 
CROS  by  FPSPEC  as  we  have  already  described  in  Section  II.  We  had  anticipated 
difficulty  in  fitting  the  gamma  spectrum  because  it  exhibits  growth  followed 
by  decay,  as  illustrated  in  fig.  7.  This  specific  problem  did  not  materialize, 
but  on  the  other  hand  we  had  much  more  difficulty  in  all  of  the  fits  than  we 
had  anticipated. 

* 2 

The  fits  are  all  weighted  so  as  to  define  x In  terms  of  the  relative 
error  at  each  decay  time.  We  ignore  the  overwhelming  covariance  introduced 
by  the  CINDER  calculation  and  use  a diagonal  weight  matrix.  We  find  that  the 
design  matrix  is  linear  with  respect  to  the  amplitudes  agi  but  nonlinear  with 
respect  to  the  decay  constants  X , . 


Mathematical  details  of  the  fitting  process  are  summarized  in  appendix  A. 
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The  nonlinear  nature  of  the  fit  isplies  that  there  may  be  more  than  one 
relative  minimum  inx^»  which  we  find  to  be  a pervasive  problem  for  more  than 
three  terms  in  the  fit.  Because  of  the  large  number  of  groups  to  be  fitted* 
we  made  the  search  for  absolute  minima  a largely  automatic  process  that  is 
built  into  FPSPFT.  We  shall  outline  the  automatic  procedure  here  and  discuss 
it  in  more  detail  in  appendix  A.  As  we  shall  see  in  Section  VI. 3,  however, 
the  automatic  search  proved  inadequate  for  fitting  the  photon  spectra,  so  that 
we  were  forced  to  adopt  additional  nonautomatic  measures  to  obtain  acceptable 
fits. 

The  automatic  search  in  FPSPFT  is  conducted  on  two  levels.  Every  fit 
begins  with  a trial  set  of  decay  constants,  from  which  a straightforward  cal- 
culation gives  the  corresponding  amplitudes  for  an  absolute  minimum  in 
The  complete  parameter  set  is  then  iterated,  following  the  multidimensional 
gradient  of  x^»  until  it  either  converges  to  a relative  minimum  or  else  clear- 
ly falls  to  converge  at  all.  The  central  problem,  then,  is  to  find  a trial 
set  of  decay  constants  that  lies  dose  enough  to  the  optimal  set  for  the  ini- 
tial gradient  to  guide  the  solution  to  it. 

In  order  to  explore  parameter  space  reasonably  thoroughly,  FPSPFT  uses 
up  to  six  different  trial  sets.  The  fitting  algorithm  always  starts  with  the 
lowest-energy  group.  Once  a group  has  been  fitted,  its  parameters  are 
used  as  the  first  guess  for  the  next  higher  group.  This  guess  is  almost  always 
successful  for  the  smooth  electron  spectrum,  at  least  in  the  sense  that  it 
finds  a corresponding  relative  minimum  (whether  or  not  that  proves  to  be  an 
absolute  minimum) . It  frequently  falls  for  groups  near  the  high-energy  end 
of  the  electron  spectrum,  however,  and  usually  falls  for  the  much  more  irreg- 
ular photon  spectrum. 

The  second  obvious  estimate  is  found  by  making  an  empirical  search  for 
regions  of  the  decay  curve  that  appear  to  be  dominated  by  a single  decay  con- 
stant. At  late  decay  times  in  the  highest- energy  groups  the  decay  is  fre- 
quently dominated  by  a single  isotope,  so  this  procedure  is  often  fruitful. 

For  the  dense  data  sets  the  search  is  made  by  fitting  a five-term  Lagrange 
interpolating  polynomial  to  the  logarithm  of  the  group  Intensity  and  examin- 
ing the  derivatives  of  that  polynomial.  For  the  sparse  data  sets  this  fit  is 
too  erratic  to  be  useful,  and  the  empirical  search  cannot  be  used. 

Four  other  trial  sets  are  based  on  the  apparent  decay  constants  at  the 
beginning  and  end  of  the  time  interval  that  Is  being  fitted.  For  the  dense 
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sets  these  are  taken  from  the  Lagrange  fits.  For  the  sparse  sets  the  loga- 
rithmic slopes  between  the  first  two  and  the  last  two  points  are  used  instead, 

* 

with  an  optional  procedure  to  extrapolate  the  estimate  averaged  over  the  last 
interval  to  the  latest  time  in  that  interval.  The  four  trial  sets  are  then 
formed  by  using  various  systematic  subdivisions  of  the  interval  between  these 
two  extreme  slopes.  In  appendix  A these  are  referred  to  as  equispaced  sub- 
divisions. 

FPSPFT  also  offers  the  option  of  dictating  initial  estimates  of  the  decay 
constants  from  input  cards,  for  selected  groups  or  for  all  groups.  In  this 
event  an  additional  option  allows  the  minimization  of  to  be  carried  out 
solely  by  choosing  the  best  amplitudes,  or  alternatively  by  iterating  on  the 
entire  set  of  parameters . 

2.  PROPERTIES  OF  INDIVIDUAL  FITS 

Clearly  the  effect  of  increasing  the  number  of  terms  m in  the  fit  is  to 
allow  the  fit  to  follow  more  of  the  undulations  in  the  decay  curves.  Figure 
12  illustrates  the  resulting  improvement  for  the  electron  spectrum  in  case  49F. 
The  bottom  portion  of  the  figure  refers  to  the  same  group  of  electrons  that 
was  displayed  in  figs.  8 and  10.  We  see  that  the  fit  crosses  the  data  at 
precisely  2m  points,  as  expected,  and  also  that  the  error  nowhere  exceeds  5% 
for  m > 3.  If  we  sum  the  fit  over  all  groups  and  compare  the  sum  to  the  sum 
of  the  input  data,  we  obtain  the  results  in  the  top  portion  of  the  figure.  We 
observe  that  the  errors  in  the  total  yield  are  roughly  equal  to  those  in  the 
individual  group  yields;  that  is,  that  no  major  reduction  results  from  averag- 
ing over  many  groups. 

The  reason  for  the  latter  result  should  be  clear  from  fig.  13,  which  shows 
the  error  at  10  seconds  after  fission  for  each  group  in  the  same  electron 
spectrum.  There  is  enough  similarity  in  the  shape  of  the  decay  curve  for 
each  group  that  many  adjacent  groups  will  have  errors  of  the  same  sign  and 
comparable  magnitude  at  any  given  decay  time.  Naturally  the  locations  of  the 
regions  of  large  error  depend  on  the  number  of  terms  fitted.  Figure  14  shows 
the  corresponding  errors  for  the  photon  spectrum.  Here  the  errors  in  adjacent 
groups  are  less  strongly  correlated,  but  the  magnitudes  display  greater  extremes. 
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In  Section  VI. 3 we  shall  refer  to  this  as  the 


L option. 
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Figure  13.  Error 
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Figure  14.  Error  in  the  Fit  to  the  17-Group  Photon  Spectrum  10  Seconds 
after  Fission  of  239pu  induced  by  Fission  Neutrons > as  a 
Function  of  the  Humber  of  Terms  in  the  Fit  (m; . The  error  for 
6 terms  is  too  email  to  plot  here. 
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3.  ADDITIONAL  MEASURES  TO  FIND  BEST  FITS 
The  measures  that  we  described  in  Section  VI. 1 were  developed  using  the 

| electron  spectrum  as  a test  problem.  For  each  of  our  six  cases  we  punched 

output  cards  containing  the  best  parameter  sets,  and  their  root-mean-square 
(ms)  errors,  for  values  of  m from  2 through  5 (2  through  6 for  the  dense 
sets) . In  a few  instances  we  had  to  exercize  the  L option  mentioned  in  Section 
VI. 1 in  order  to  get  the  best  fit. 

2 

In  an  attempt  to  broaden  the  search  for  absolute  minima  in  x . we  then 

used  each  of  these  punched  parameter  sets  as  the  initial  estimate  for  the  other 

five  energy-isotope  cases,  allowing  FPSPFT  to  iterate  all  parameters  from  chat 

2 

start.  In  none  of  these  cases  did  we  find  a further  improvement  in  X » which 
we  found  to  vary  smoothly  from  group  to  group,  so  we  accept  these  parameters 
as  giving  the  best  possible  fits. 

To  our  dismay,  in  both  the  17-  and  38-group  photon  spectra  we  found  that 
the  L option  gave  better  fits  in  some  groups  and  worse  fits  ip  others.  Ac- 
cordingly, we  had  to  include  Doth  alternatives  in  the  cross-fertilization 

process  described  above.  This  t)«ae  the  cross-tertilization  produced  many 
2 

improvements  ir.  x » some  of  theA  very  large  improvements,  especially  in  the 
4-  and  5-tern  fits.  Some  of  these  triggered  improvements  in  adjacent  groups 
also.  When  the  improved  solutions  were  recycled  into  the  crcss-fertilization 
they  produced  still  further  improvements.  This  process  finally  converged 
after  four  iterations.  It  required  much  clumsy  hand  work  with  thousands  of 
cards,  so  that  we  were  eventually  driven  to  writing  a separate  computer  prog- 
ram just  to  read  and  compare  the  rms-error  cards. 

4.  OVERALL  FITTING  ERRORS 

For  each  case  FPSPFT  calculates  rms  relative  errors  routinely  for  the  fit 

to  each  group  in  a spectrum,  for  the  ensemble  of  all  points  in  all  groups,  and 

for  the  summed  fits  compared  to  the  CINDER/FPSPEC  sums.  In  support  of  our 
qualitative  remarks  that  are  illustrated  by  figs.  13  and  14,  the  rms  error  in 
the  suos  is  always  smaller  than  that  for  the  ensemble,  but  usually  not  very 

much  smaller.  For  practical  purposes  in  judging  the  effect  of  using  a partic- 

ular number  of  terms  in  the  fit,  it  is  more  useful  to  calculate  a yisld- 

2 

weighted  rmo  error,  defined  as  the  cquare  root  of  the  weighted  average  of  x > 

3 

vhere  each  group  is  weighted  by  its  integrated  yield  from  C to  10  seconds. 

In  addition,  m compounding  these  results  over  our  nix  energy-isotope  cases 
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it  is  appropriate  to  weight  the  dense  and  sparse  cases  by  the  number  of  times 
* 

in  the  case. 

The  resulting  point-and-yield-weighted  rms  error  (PYWE)  provides  a con- 
venient index  for  evaluating  the  penalty  associated  with  using  a given  number 
of  terms  in  the  fit.  The  result  is  plotted  as  the  solid  lines  in  fig.  15.  We 
see  that  for  only  two  terms  the  PYWE  is  about  30%  for  both  electrons  and  pho- 
tons regardless  of  the  number  of  groups  in  the  spectrum.  The  PYWE  is  less 
than  2%  for  5 terms  and  less  than  1%  for  the  three  cases  in  which  6 terms  can 
be  fitted,  although  in  the  latter  cases  there  is  more  than  a factor  of  2 
spread  between  the  20-group  beta  and  38-group  gamma  PYWEs.  Figure  15  in  con- 
junction with  the  error  estimates  for  the  ENDF/B  data  base  discussed  in  Sec- 
tion IV  thus  allows  an  informed  choice  of  how  many  terms  to  use  for  a partic- 
ular application. 

5.  EXTRAPOLATION 

Fitting  a decay  curve  with  an  exponential  series  permits  extrapolation 
to  both  earlier  and  later  times.  Clearly,  the  fitting  process  cannot  usually 
pick  up  components  with  decay  constants  whose  reciprocals  are  much  smaller 
than  the  earliest  or  much  larger  than  the  latest  time  included  In  the  fit. 

There  are  so  few  beta-decay  half-lives  less  than  0.1  second  that  we  can 
be  certain  that  CINDER  calculations  will  extrapolate  reliably  to  zero  time; 

On  the  basis  of  the  half-lives  in  ENDF/B  we  find  that  less  than  1%  of  the 
integrated  electron  yield  from  zero  to  infinity  occurs  before  0.1  second.  On 
the  other  hand,  we  have  pointed  out  in  Section  I that  we  have  deliberately 
Ignored  the  gamma  emission  at  early  times.  This  emission,  which  accounts  for 
roughly  45%  of  the  photon  emission  from  zero  to  infinite  time,  includes  a very 
iatense  prompt  burst  within  the  first  nanosecond  followed  by  decay  of  more 
than  80  known  activities  with  half  lives  between  a few  ns  and  80  Ms  (ref.  13). 
Tne  delayed  component,  which  is  usually  attributed  to  isomers  of  the  primary 
fission  fr.igrr.en c a , becomes  negligible  relative  to  the  photons  following  b ec» 
decay  at  a decay  tin*-  of  •’bout  s.  millisecond. 


Five  terms  with  two  parameters  per  term  permit  an  fit  to  the  &~p°int 

sparse  sets  in  most  groups,  but  the  fits  are  physically  unrealistic.  Hence 
we  shoulu  reduce  the  weight  of  the  sparse  sets. 
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Since  gamma  emission  more  than  a thousand  seconds  after  fission  is  almost 
entirely  from  de-excitation  of  the  daughters  of  beta  decay,  both  emissions 
should  have  comparable  inaccuracies  in  extrapolating  the  total  emission  to 

later  times.  To  illustrate  such  extrapolation,  in  case  49F  our  fit  from  0.1 

3 

to  10  seconds  predicts  that  84%  of  the  total  number  of  electrons  that  will 

3 

be  emitted  have  emerged  by  10  seconds,  and  by  extrapolation  that  100%  will 
4 

have  emerged  by  10  seconds.  Unfortunately,  if  we  extend  the  fitting  range 

t ^ 

to  10  seconds  we  find  that  only  70%  of  the  emission  is  actually  complete  by 

3 4 4 

10  seconds,  13%  more  emerges  by  10  seconds,  and  extrapolation  beyond  10 

seconds  predicts  at  least  an  additional  5%  to  infinity. 

In  fig.  5 we  showed  selected  20-group  electron  spectra  for  case  25F  out 

4 3 

to  10  seconds  after  fission.  If  the  fit  stops  at  10  seconds,  no  group  ex- 

4 

trapolates  to  an  intensity  at  10  seconds  that  is  greater  than  the  bottom  of 
the  figure  [10  electrons/ (MeV*s* fission) ] , although  all  groups  below  5 MeV 

should  exceed  that  limit.  In  fig.  6,  on  the  other  hand,  we  see  11  out  of  38 

3 4 

photon  groups  fitted  to  only  10  seconds  that  extrapolate  to  values  at  10 

seconds  that  are  greater  than  the  bottom  of  the  plot  [10  pho tons /(MeV’S' fis- 

sion) ] , including  one  that  extrapolates  to  a greater  intensity  than  the  cor- 
rect value  from  CINDER. 

To  summarize,  roughly  30%  of  the  beta-decay-dominauud  emission  occurs 
3 3 

after  10  seconds.  The  intensities  at  10  seconds  and  integrals  up  tc  that 

3 

time  can  be  represented  accurately  by  parameters  fitted  from  0.1  to  10'  seconds, 
but  extrapolating  that  fit  misses  about  5%  of  the  total  yield  to  infinity. 

The  extrapolated  intensity  is  almost  invariably  too  small  by  at  least  a factor 

4 

of  ten  at  10  seconds,  and  the  error  is  much  more  erratic  in  the  photon  spec- 
trum than  in  the  electron  spectrum. 

6.  JOINT  FITS 

2 

The  existence  of  multiple  relative  minima  in  x suggests  that  each  group 
in  a spectrum  might  be  fitted  reasonably  accurately  by  a vide  variety  of  decay 
constants,  as  long  as  the  amplitudes  are  chosen  accordingly.  This  suggests 
in  turn  that  a single  set  of  decay  constants  might  be  used  for  all  six  energy- 
isotope  cases  simultaneously,  thus  reducing  by  nearly  a factor  of  two  the 
total  number  of  parameters  required  to  represent  all  of  the  data  with  a spec- 
ified accuracy.  If  such  a joint  fit  is  successful,  it  offers  an  important 
but  less  obvious  additional  advantage.  The  sparse  data  sets  are  divided  very 
unequally,  with  time  intervals  of  a factor  of  five  alternating  with  a factor 


of  two,  We  find  that  the  overdetermined  5-tero  fits  to  these  sets  oscillate 
systematically  towards  and  away  from  the  fits  to  the  dense  sets,  taking  advan- 
tage of  the  factor-of-five  gap  in  which  they  are  unconstrained  to  make  it 
easy  to  fit  the  pairs  of  points  that  are  only  a factor  of  two  apart.  Thus, 
the  use  of  a single  set  of  decay  constants  (derived  in  this  case  only  from  the 
dense  sets)  probably  insures  a fit  that  is  closer  to  the  true  (but  uncalcu- 
lated) values  in  the  long  gaps. 

We  did  not  recognize  the  desirability  of  combined  fits  soon  enough  to 

design  simultaneous  fitting  of  many  cases  into  FPSPFT.  Accordingly,  we  tried 

* 

the  simple  alternative  of  averaging  together  the  corresponding  decay  constants 
from  the  six  cases.  Although  the  rms  difference  of  the  six  individual  values 
from  their  average  is  often  of  the  order  of  many  tens  of  percent,  the  joint 
fits  display  an  unexpectedly  small  Increase  in  the  PYWE  defined  in  Section 
VI. 4.  In  an  attempt  to  explore  this  further,  for  both  photon  group  struc- 
tures we  added  an  additional  step  in  finding  the  joint  fits.  We  used  the 
average  decay  constants  from  the  free  fits  as  initial  estimates  for  each  case 
in  turn,  but  allowed  FPSPFT  to  iterate  both  amplitudes  and  decay  constants 
away  from  those  estimates.  We  then  re-averaged  the  resulting  decay  constants 
for  the  final  parameter  set.  The  resulting  Improvement  was  small.  In  a few 
groups  the  composite  fit  was  actually  worse  than  before,  so  that  we  had  to 
return  to  the  original  average.  # 

The  dashed  lines  in  fig.  15  show  the  values  of  PYWE  for  the  final  joint 
fits  compared  to  the  separate  (and  therefore  best)  fits  given  by  the  solid 
lines.  In  examining  this  figure,  we  find  empirically  that  the  PYWE  of  the 
joint  fit  is  approximately  0.6%  greater  than  for  the  best  fit  to  the  20-group 
electron  spectrum,  roughly  Independent  of  the  number  of  terms  in  the  fit.  The 
corresponding  differences  for  the  17-  and  38-group  photon  spectra  are  1.3%  and 
2.0%,  respectively.  We  note  also  that  in  many  cases  the  joint  fit  with  m 
terms  has  both  a smaller  PYWE  and  fewer  parameters  than  does  the  best  fit 
with  m-1  terms,  in  addition  to  giving  a physically  superior  fit  to  the  sparse 
data  sets. 


*As  noted  in  appendix  A,  FPSPFT  stores  the  decay  constants  in  order  of 
decreasing  size.  Corresponding  constants  are  those  with  the  same  index. 
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SECTION  VII 

SUGGESTIONS  FOR  FURTHER  WORK 


It  is  clear  from  fig.  11  that  the  gradual  accumulation  of  a more  complete 
data  base  has  made  a very  large  difference  in  fission-product  calculations 
since  Dieckhoner’s  work  was  completed.  We  have  also  emphasized  in  Section 
III  that  on  the  average  less  than  half  of  the  total  energy  emission  at  times 
less  than  one  minute  after  fission  is  accounted  for  by  detailed  spectra  in 
ENDF/B-4,  and  in  Section  V.5  that  fits  should  extend  to  at  least  10^  seconds 
if  reliable  integrals  to  infinity  are  needed. 

Version  5 of  ENDF/B  is  scheduled  to  have  more  than  1,000  fission-product 
nuclides  in  its  decay-property  tables,  including  about  240  for  which  detailed 
spectral  information  will  be  included.  As  we  pointed  out  in  Section  II,  the 
new  version  will  also  include  nuclldlc  yields  of  fission  products  that  have 
been  calculated  using  improved  models  for  the  pairing  effect  and  the  branching 
fraction  between  ground  and  isomeric  states.  Thus,  it  would  be  attractive 
to  repeat  the  fitting  of  electron  and  photon  spectra  after  Version  S has 
finished  its  lengthy  review  and  testing. 

The  computational  tools  for  calculating  the  spectra  are  already  complete. 
The  fitting  code  FPSPFT  can  already  handle  extension  beyond  10"*  seconds,  but 
needs  to  have  an  additional  algorithm  added  that  would  perturb  solutions  after 
convergence  in  a systematic  way  in  an  extended  effort  to  locate  all  of  the 
local  mimima  in  x • If  experience  by  then  has  shown  that  AFWL’s  particular 
applications  benefit  from  the  more  compact  parameter  set  that  results  from 
using  a common  set  of  decay  constants  for  all  isotope-energy  cases,  it  would 
be  desirable  to  add  rigorous  simultaneous  fitting  of  all  cases  to  FPSPFT. 

The  spectrum  of  photons  from  the  prompt  (<  1 ns)  burst  has  been  measured 

for  a number  of  cases,  but  the  subsequent  decay  (1  ns  to  1 ms)  of  lsoeterlc 

states  has  been  measured  carefullv  only  for  well-moderated  neutrons  Incident 
235  239 

on  U and  Pu.  If  this  fragmentary  information  Is  useful,  parameter isat ion 
in  any  given  group  structure  would  be  relatively,  straightforward. 
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APPENDIX  A 

MATHEMATICAL  DETAILS  OF  FITTING  DECAY  CURVES 


? 

| 

\ 


In  this  appendix  we  shall  designate  all  vectors  by  unsubscripted  lower- 
case letters,  all  matrices  by  capital  letters,  and  all  scalars  by  lower-case 
Greek  letters.  Our  treatment  basically  follows  that  of  Hamilton  (ref.  19). 

For  each  group  in  a given  spectrum,  CINDER  and  FPSPEC  give  us  a table  of 
intensities  y^  as  a function  of  decay  times  t^.  For  convenience,  we  treat 
these  as  vectors  y and  t,  respectively.  The  covariances  <dy^  dy^>  of  the  y^ 

are  embodied  in  a covariance  matrix  V.  As  mentioned  In  Section  VI. 1,  we  have 

2 

taken  V to  be  diagonal  with  elements  V^  * (0.01  y^)  ; that  is,  we  have  ignored 

correlations  induced  by  the  fact  that  our  "data"  are  calculated  from  a common 

data  base  (as  opposed  to  being  independent  measurements) , and  have  arbitrarily 

assumed  a 1%  standard  deviation.  The  weight  matrix  W = V ^ then  causes  the 

2 

fit  to  be  weighted  by  the  relative  error,  with  the  generalized  x in  units  of 
2 

percent  . 

Our  task  is  to  find  m decay  constants  X ^ and  an  equal  number  of  amplitudes 
(Xj  such  that  the  values  of  y^,  defined  by  the  relation 


m 


-Z 

>1 


a. 


-'iS 


(Al) 


give  a best  fit,  in  the  least-squares  sense, 
convenience  we  shall  collect  the  parameters  a 
sional  parameter  vector  p such  that 


to  the  FPSPEC  values 


For 


j and  X into  a single  2m-dimen- 


PJ  "aJ 


' j-m 


0 < j - m 
m < j - 2m 


We  put  the  amplitudes  first  in  p for  convenience  in  a preliminary  linear  fit 
using  ass'-  ed  initial  values  for  the  decay  constants,  as  discussed  below.  All 
infwrmation  concerning  the  times  and  the  structure  of  equation  Al  is  then 
embodied  in  the  design  matrix  D with  elements 
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D«  5 


_tiXj 

a +• 


0 < j = m 


(A2a) 


-t,  a.  e 
i j-m 


-t  A . 
i j-m 


m < j - 2m  . 


(A2b) 


In  the  usual  fashion  for  covariant  fitting,  the  quantity  to  be  minimized  is 

X2  i d • Wd  - 3 W d , 

where  d = y - y is  the  difference  between  the  fit  and  the  data,  and  3 denotes 

2 

the  transpose  of  d.  We  minimize  x by  finding  the  value  of  the  parameter 

2 

vector  p for  which  the  gradient  of  x with  respect  to  p vanishes: 

Vx2  -25  Wd  - 2 5W  (y-y)  - 0 , 

and  hence  6 W y » D W y . (A3) 

2 

Here  we  have  used  the  chain  rule  to  differentiate  x with  respect  to  the  ele- 
ments p^.  Equation  A3  is  the  basic  form  of  what  is  traditionally  called  the 
normal  equation. 

Up  to  this  point  the  solution  is  completely  general  and  we  have  ignored 
the  actual  form  of  the  design  matrix  in  equations  A2a  and  A2b.  if  the  elements 
of  the  design  matrix  are  independent  of  the  fitted  parameters,  equation  A3  re- 
duces to  a set  of  linear  equations  with  an  unique  solution.  Since  this  is  not 
the  case  for  equation  A2,  the  problem  is  nonlinear  and  must  be  solved  by  iter- 
ation. We  note,  however,  that  the  amplitudes  a.  do  not  appear  in  equation 

A2a,  so  that  if  we  assume  trial  values  of  the  \ . we  can  immediately  deduce  the 

2 J 

amplitudes  that  will  minimize  x for  that  particular  set  of  decay  constants, 
using  only  the  parts  of  p and  D that  correspond  to  a.  Thus,  we  can  begin  the 
solution  by  estimating  only  the  decay  constants 

If  we  designate  this  initial  set  of  trial  parameters  by  p and  the  cor- 
responding initial  fit  by  y,  the  next  step  is  to  seek  an  improved  parameter 
set  p of  the  form 

p - p + q , (A4) 

2 

in  which  q will  be  chosen  to  minimize  x • This  produces  a corresponding 
change  in  the  fit  given  by 


y s y + D q , 


(A5) 


where  we  have  kept  only  the  first  terms  in  a MacLaurin  series  in 
change  q,  and  taken  advantage  of  the  fact  that  the  first  partial 
of  y with  respect  to  the  parameters  are  just  the  elements  of  D. 
this  approximation  we  have  linearized  the  normal  equation,  which 

6 W (y  fDq)  - 5 W y , 

and  hence  (CWD)  q ■ B W (y-y) 

The  product  N I B W D is  known  as  the  normal  matrix.  From  its 
can  readily  find  a solution  for  q, 

q • N'1  5 W (y-y)  , (A6) 

in  which  the  difference  y - y is  simply  the  error  in  the  trial  fit.  Equation 

A6  is  not  rigorous  because  expansion  A5  is  not  exact.  However,  if  q is  small 

compared  to  p we  can  expect  to  find  that  the  new  value  of  p calculated  from 

2 

equation  A4  produces  a reduction  in  x , so  that  by  recalculating  D using  the 

new  parameters  we  can  Iterate  until  q becomes  vanishingly  small. 

2 

If  the  first  iteration  Increases  V we  know  that  the  initial  estimate 
for  p was  too  far  in  error  for  equation  A5  to  be  a satisfactory  approximation. 
Accordingly,  in  the  hope  that  the  multidimensional  direction  of  q is  more 
nearly  correct  than  its  length,  we  can  replace  equation  A4  with 

P - p + q 6 , (A7) 

where  6 is  a damping  constant.  FPSPFT  tries  values  of  6 beginning  at  1 and 

2 -*5 

decreasing  by  a factor  of  2 until  a decrease  in  x occurs  or  until  6 < 10  . 

If  an  Improvement  does  occur,  FPSPFT  recalculates  D and  q from  the  successful 

parameters,  doubles  the  previous  value  of  6,  and  tries  again.  Usually  it  can 

continue  doubling  in  each  subsequent  iteration  until  the  value  of  6 returns  to 

unity. 

When  q is  large,  it  may  produce  such  large  changes  in  p that  one  of  the 

decay  constants  will  be  displaced  past  one  of  its  neighbors.  If  this  happens, 

2 

even  if  a decrease  in  x would  result,  it  Is  treated  as  an  unsuccessful  iter- 
ation,  and  the  damping  constant  is  reduced  successively  until  each  decay  con- 
stant is  at  least  10%  greater  than  the  one  below  it  (reduced  to  5X  after  five 
iterations) • 


the  small 
derivatives 
By  making 
now  reads 


inverse  we 
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FPSPFT  stops  at  a preset  number  of  iterations  (usually  25)  unless  the 

solution  converges  sooner.  The  criterion  for  convergence  is  that  x decrease 

by  less  than  a small  fraction  £ of  its  previous  value.  After  the  first  iter- 

—8 

ation  (for  which  any  decrease  is  acceptable),  £ is  set  to  10  . Thereafter 

it  Is  doubled  in  each  Iteration,  so  that  the  definition  of  convergence  becomes 

h 

increasingly  less  restrictive.  This  procedure  allows  FPSPFT  to  work  with 

slightly  overdetermined  cases  without  wasting  an  inordinate  amount  of  time 
2 

trying  to  get  x ever  closer  to  fcero.  It  also  terminates  a fruitless  search 
for  a minimum  that  is  too  poorly  defined  to  merit  accurate  parameters. 

In  Section  VI. 1 we  discussed  two  methods  of  generating  initial  estimates 
of  the  decay  parameters  that  are  required  to  begin  the  iterative  search.  One 
of  these  is  to  use  the  solution  from  the  previous  group,  and  the  other  makes 
use  of  the  empirical  logarithmic  slope  of  the  decay  curve.  The  latter  method, 
which  we  attempt  only  for  the  dense  data  sets,  requires  an  automatic  means  of 
defining  regions  in  the  decay  curve  for  which  the  logarithmic  slope 
A = -d(log  y)/dt  of  the  intensity  y is  nearly  constant  for  an  extended  in- 
terval. This  means  that  we  seek  regions  in  which  the  second  and  third  deri- 
vatives of  log  y are  both  small  compared  to  log  y Itself.  We  have  found  by 
trial  and  error  that  the  product  of  the  normalized  second  and  third  deriva- 
tives is  a simple  and  satisfactory  measure  of  this  property,  so  FPSPFT  simply 
chooses  values  of  A that  correspond  to  minima  in  this  product,  beginning  with 
the  least  positive  minimum.  The  selected  values  are  then  sorted  into  descend- 
ing order. 

The  four  remaining  methods  that  FPSPFT  uses  to  find  Initial  estimates 
produce  equispaced  values  in  a relatively  blind  attempt  to  explore  the  pa- 
rameter space  within  which  the  decay  constants  are  expected  to  lie.  We  begin 
by  estimating  the  logarithmic  slopes  at  the  beginning  and  end  of  the  time 

3 

interval  0.1  to  10  seconds,  which  we  term  A^  and  A#,  respectively.  For  the 
dense  sets  we  take  these  directly  from  the  empirical  analysis.  For  the  sparse 
sets  we  cake 

Xfa  - logCyj/y^/Ctj-tp 


*Sinca  225  x 10“ 8 - o.34,  there  la  clearly  no  reaaon  to  attempt  more  than  25 
iterations  under  this  definition  of  convergence. 
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f 

i 

\ 

i 

l 


and 


X 

e 


(A8) 


where  there  are  n decay-times  in  the  data  set.  FPSPFT  includes  an  option 
to  extrapolate  X g approximately  to  tR,  rather  than  using  the  average  of  tn  ^ 
and  tn  implied  by  equation  A8.  Setting  r = Xg/Xb,  we  then  define  three  con- 
stants: 


6 

w 

6 

m 


rl/(m-l) 
rl/  m 


wide 

medium 


rl/  (m-fl)  , narrow 

where  m is  the  number  of  terms  in  the  fit  and  the  subscripts  stand  for  the 
indicated  adjectives.  These  constants  are  used  to  generate  the  four  geometric 
series  of  decay  constants  shown  in  table  Al.  The  wide  equispaced  set  spans 


» 


Table  A1 

EQUISPACED  ESTIMATES  OF  DECAY  CONSTANTS 

£?■  utn< 


Index 

High 

Wide 

Na  cror* 

Low 

1 

X, 

b m 

X,  6 
h n 

Xb 

2 

X,  o v 
b m 

X.  6 
b w 

X.  6 2 
b n 

Xb  5m 

3 

X.  f 3 
b a 

X,  o 
b w 

X.  6 5 

b tt 

X.  6 2 
b m 

« 

• 

• 

• 

• 

* 

• 

• 

* 

• 

• 

o-l 

X /6 
e m 

• 

X /5 
c w 

• 

X /6  2 

s n 

« 

X /<$  2 
e in 

m 

X 

e 

\ 

v» 

X /<S 
e n 

X /6 
e m 

In  Section  VI. 3 we  refer  to  this  as  the  l Option, 
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the  interval  from  X^  to  in  equal  logarithmic  steps.  The  narrow  equispaced 
set  is  compressed  so  that  it  is  wholly  interior  to  the  same  interval.  The 
high  and  low  equispaced  sets  u3e  an  intermediate  spacing,  and  are  offset  so  as 
to  end  on  Ag  or  begin  on  X^,  respectively.  FPSPFT  does  not  use  estimates  that 
lie  outside  the  empirical  range,  since  we  have  found  that  it  is  easy  for  the 
iteration  process  to  move  parameters  outwards  from  within  the  range  but  some- 
times is  difficult  to  move  inwards  from  outside  the  range. 
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APPENDIX  B 

FORMAT  OF  TRANSMITTAL  CARDS 

We  have  used  two  different  formats  for  the  ouneheu  cards  that  contain 
the  parameter  set3  Chat  describe  the  fitted  decay  curves.  The  first  form  is 
used  for  the  free  fits  (or  best  fita),  which  have  different  decay  constants 
for  each  energy-isotone  case.  The  aata  are  organized  into  self-contained 
decks,  of  which  an  example  is  shown  in  fig.  Bl  for  the  6-tera  fit  to  the  17- 
group  photon  spectrum  from  case  49H.  Most  of  the  format  is  explained  by  an- 
notations in  the  figure.  The  information  at  the  end  of  the  deck-heading  card 
correlates  the  carda  with  their  corresponding  printed  output.  The  primary 
information  is  a succession  of  pairs  of  parameters  (A^cep  for  the  successive 
terms  In  equation  Al.  A separate  sub tab  'o  in  this  form  is  given  for  each 
group . 

The  second  form  of  transmittal  cards  is  used  for  the  joint  fits,  in  which 
a single  set  of  decay  constants  serves  all  six  energy-isotope  cases.  Since 
this  parameter  set  emphasizes  compactness,  we  have  used  a card  format  that 
minimizes  wasted  space  on  the  carda.  Figure  B2  shows  the  beginning  and  end 
of  the  17-group  photon  parameter  set,  together  with  an  intermediate  portion 
of  the  set.  In  this  format  we  have  used  descriptive  heeding  cards  to  make  the 
format  as  nearly  self-explanatory  as  possible.  Here  the  boundaries  of  the 
group  structure,  which  is  constant  for  the  entire  set,  are  given  first.  Next 
we  give  the  decay  cons  cants  for  all  2-term  fits,  followed  by  separate  sub- 
tahles  of  the  amplitudes  for  each  2-tei.'m  fit,  and  so  on  for  each  number  of 
terms  in  succession.  Consequently,  the  primary  information  consists  of  a 
list  of  values  (not  pairs  of  values)  of  a particular  type  of  parameter  for 
successive  groups  in  the  spectrum. 

In  both  forms  all  of  the  significant  integer  information  is  given  in 
fields  with  Fortran  format  15  that  end  on  a multiple  of  5 columns.  The  energy 
widths  and  boundaries  of  the  group  structure  are  in  F5.2  format  in  fields  that 
also  end  on  multiples  of  5 columns.  All  other  floating-point  values  are  in 
E10.3  fields  that  end  on  a multiple  of  10  columns. 
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9U 

1 


239 

<>,25 


Atomic  number  of  fissioning  nucleus 
— Mass  n « it  it 

Neutron  energy  ( 1 - thermal,  2 - fission,  3 • 1U  MeV  ) 

- Type  of  particle  ( 1 • electron,  2 ■ photon  ) 

Number  of  groups  in  spectrum 

r— — — — — Number  of  terms  fitted 

Y ^ ■ Date  and  time  cards  were  punched 

LAMBDA 


[ 


3 

.75 


17  6 6/12/77 

Group  1 upper  energy  and  width  (MeV) 


3.108E-00  l,049E-««  1 , 31  IE*00  5.076F-&5  1.813F-0I  8.969E-05  b,08i£.02  3,lS9E-<ib 

V » a,  — i 


au 


1.922E-02  3.710E-06  l,2bbfc-02  3.  177E-06  V— *3  *—\h 

2 5.50  ,50 

3,229E>00  9,755E*05  i,233F*0fl  4.203E-05  1.788E-01  7.801E-05  J,011E*02  2,b9lE-0b 
l,l«5E«02  b, 36bE-0b  4.034E-0S  l,810E-0o 

3 5,00  ,50-^ Group  3 upper  energy  and  width 

S,210E*00  7,113E-«5  l,2bbE*00  3,2«bE-05  1.770E-01  b,HiaE-05  1 ,ttf97t-02*i,507E-«ib 
9.901E-03  7.278E-0S  3.899E-03  2,l«bE*0b 
a 4,50  ,50 

2, 8 1 2E*00  2.124E-04  7.886E-01  5,b?0£-05  4.782E-02  5,75bE-05  1.471E-02  b,9bbE-05 
7.782E-03  2,465E«05  3.598E-03  1.202E-05 

5 4 00  50 

3,082Et'J0  2l«95E-0«  i,098C*00  l,23uf-04  3,929F.-02  1,149E*04  i.375fe-02  l,2b3E-H4 
5.244E-03  1,720E*05  t,t30E-03  3,46bE-0b 

6 3,50  ,50 

2,930Ef00  2.392E-04  9.533E-01  7,b23E-05  6.398E-02  7.312E-05  1.749E-02  b,997E-05 
6, 38bE-0  3 4.597E-05  1,561E*03  9,690E-0b 

7 3,00  ,50 

2, 369Et00  l,284E«03  3, 1 8bE-0 1 • 1 , 763E-04  b,4!0E-02  7,3b8t-04  1.B90E-02  3,323E*0« 
8,00bE-03  1.152E-04  9,b99E-04  1.382E-05 

8 2,50  ,50 

3,207E*00  9.103E-04  1,225E*00  4,552E-04  9,55lE-02  3.375E-04  2.791E-02  3,«87t-04 
9,b0bE-03  2.127E-04  7.210E-04  2.924E-05 

9 2 00  50 

2,925E>00  J*559f -03  8,«93E*0l  l,07«E-03  1.052E-01  1.408E-03  2.731E-02  1.M43E-03 
9.320E-03  3.190E-04  6,7871-04  5.263E-05 
10  1,50  ,50 

2.815E+00  1 , b2 IE-02  7,0b7E*0i  4,3«9E-03  l,l4bE-01  b,835E*03  2,889E-02  2,83lE-03 
9.434E-03  1.573E-03  8.897E-04  l,b7uE-04 
It  1,00  ,25 

2,2 59Ef00  b, 847E-02  b,579t-0l  7,l5lE-02  I.212F-01  «,b«3t-03  2.842E-02  t,54lE-0i 
7 , 097E-03  5.521E-04  7, 9621-04  l,9b9E-04 

12  ,75  ,25 

2,671E*00  9 , 824E-02  3,0bflE-0J  b,6lbE-02  t,H57E-0l  2,0«bE-02  2,532E-02  2,778E-03 
7, 7S4E-03  l , 557E-04  9.149C-04  t,84SE*04 

13  ,50  ,25 

2, bb5£  + 00  b, 78 7E-02  3,40iE-0i  2.957E-02  1.098E-01  2.967E-02  2,bl9E-02  3,0b3E-0J 
b, 579E-03  1 , 038E-03  6,b0bF.-0a  3.J48E-04 

14  ,25  ,10 

2,859E«00  4, 764E-03  8,052E-01  l,SUE-03  l,02bE-01  t,608E-03  2,b33E-02  l,3b0t-01 
6, 798C-03  3, 185E-04  l,2b9F-03  l,4btE-04 

15  ,15  ,85 

2, 1 76E+00  b, 363E-02  b,502E-0t  b,897E-02  9,470E-02  2,0421-03  2,8«3E-02  b,22«C-fla 
7,1581-01  7,98bE-05  7.092E-04  2.060E-05 

16  ,10  ,05 

2,365E*00  7, 209E-04  2,751E-0t-2,03«E-04  5,bt5E-02  4,035E«»4  1.486E-02  2,41*E-®« 
5,269E«03  UI26E-04  9,  163E-04  3,  14lE-®55  rtwrmw  IT  upper  emeaMTW  and  width 

1,14 1E  + 00  3, 4 12f -04  t,0«9t*00  1,821E*04  l ,06bE-0|-2, 5221-04  3,45lE»02  2,7I5E»04 


b,562E-05  1,894F*04  9,75bE-04  1.084E-04 

*-x«j  v^a* 


s 


*-*3 


Figure  Bl<  Tranamittal  Card  Format  for  Frae  Fits  (Annotated) 
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-Top  of  group  1 


17-GROuP  ENERGY  ri 
'b«25  5,50  5, a*  a, 5m 
, 05  0,00 


2.T£0h  0EC*7  COn 

1 1.9MUF-01 
3,R65t»02  3,a2Hf-02 
2, 4196-02  2,24*f-R2 

2 1.32PE-02 
1,5306-03  1 , 7 1 of -01 
2 , 36 1 E“H i t,3*tE-05 

AMPLITUDES  FOP 
t 4,6P  7f -04 

2s02at-03  4,98af-03 
1.426E-P3  5.9956-05 
2 9,9686*36 

9,323E-«S  l.ulbE-Pa 
*»,  6521*05  I,a65t-0a. 

^ Qroup  3 6 


OljNOAS  If.S 

<|,0P  3,5m  3,00  ?,5m  2,.}0  1,50 
Bottom  of  group  1? 


STANTS  F Op  E*CH  GROUP- 
1.852E-01  1,9276*01  2,3t4f-02 
4,9636-02  S.615E-0I  1.936E-01 
— Decay  constant  for  1st 
S.4I7E-03  o.PMF-03  «,«9«E.0J 
2,  1 02E“0i  3.358E-03  I,655E-«3 

92  235  Gamma  FAST 


1,00  .75  ,50 

- Applies  to  all 
17  GROUPS 
2.003E-02  5, 721E- 
l,5B9fc«0l  u,07BE* 
term  in  group  17 
2,4blt«03  3 , 78b£ • 
2,630E-03  1 , 79tJF. « 

17  GROUPS 


•,is  • ,s  ,ld  Group 
6 oases  j 


4,0086-04  3.1336*04  3.2116*04  5,7386*04  3.6RR6 
2,3696-02  2,6286-01  1,32«E-01  8.445E-02  8.425E 
-Amplitude  of  1st  tern  la  group  17  for  case 


1.3806-05  1.287E-05  7,3056-05  2,3956-05  9.177E 
0,3796-04  1.1246-03  1,3776-03  1,3026-03  2.4576 


Amplitude  (Of  2nd  term  In  group  17  far  case 


02  3,026fc«0<2  1st 

02  5,4356-01  term 

03  1.963E-03  2nd 
03  4,3296-03  tem 

2 TERMS 

• 04  2, 444)1-03  1st 
-03  2, *696-81  tem 
25F 

• 05  5 , 3 35t *05  2nd 
•04  2.056E-04  tem 

2$F 


3*TERm 

OECAV  CONSTANTS  FOR 

fACH  GROUP 

17,  GROUPS 

1st 

1 

2,9l6E*00 

2.068E-01 

4,72*6-01 

2,3!0Et00 

2,3776*00 

1 ,b2«E*00 

l , 88 1 E-0 1 

1 ,40lE«00 

l ,4«2E*0W 

5.1P9E-0! 

7,  19RE-01 

2, 3B6E-M1 

2.0766-01 

1,5806*04 

6,71  IE-01 

tarn 

2,99 16*04 

3, 1«5E*00 

l , 96P6-02 

2nd 

2 

1 , 58 |E-0 1 

1 ,4996-02 

3,29«E-02 

1.835E-02 

1, 7a3fc-02 

1.7696-02 

2.200E-02 

2.625E-P2 

2.516F-02 

3,3«lF-«2 

2, 144F-02 

2.376E-02 

3,  123E-02 

4 , 1 3 7E-02 

tem 

2,095t-«2 

6.974E-03 

l,600E-03 

3rd 

3 

1.324E-02 

4, 342E-0J 

4,5032*03 

4.307E-03 

2.336E-03 

2.627E-03 

1 , 209E-03 

1.473F-03 

! .4596-0 j 

1 , A44E-03 

1.505F-03 

l,2»7C-03 

1,5206-03 

1 .469E-03 

tem 

1 .948E-03 

4, 985E-P4 

amplitudes  for 

9?  235 

gamma  fast 

17  GROUPS  3 

TERMS 

1st 

1 

5.587E-04 

4,  181E-04 

3.915C-04 

5.347E-04 

8.083E-04 

7,21 3C-04 

2.026E-03 

3.851C-03 

6,5106-03 

3.4186-02 

3.0176-01 

1.377E-01 

9 , 080E -02 

1.I87E-02 

3,1  46t -01 

tem 

l , 6 1 8E-03 

1, J81E-04 

1.298E-03 

2nd 

2 

3,500E-0« 

1.294E-05 

3, 077E-05 

2.618F-04 

4, 778E-04 

2.O31E-04 

1.409E-01 

3.660E-03 

1,  144E-P2 

5, 31 0E-03 

5,e72E-0i 

4,87  3E-03 

6.134E-03 

5,5206-03 

tem 

1.234E-0J 

1 ,229f-04 

4 , W80E-05 

3rd 

3 

9,8136-06 

6, 9 1 5£ -06 

9,608E"-7o 

6,35 1 E-05 

2 , l 7 1 E*05 

3, 779E-05 

7 , 29 JE-05 

1 , 164E-04 

3,8666-04 

4,8506*04 

3,505F-0« 

5, 364E-04 

1.975E-0O 

3,87 1 E-05 

tem 

ft*  977t*»#5  79*bb€+*5 


amplitudes  for 

94  239 

Gamma 

1 a*MEV 

17  GROUPS  6 

TERMS 

1 

l,459g-0« 

1,31 1F-04 

7,0586-195 

2.229E-04 

3, 601E-04 

2,35 l fc-04 

1 , 33 It -0  3 

1st 

8.953C-04 

2.875E-03 

1 .661E-0? 

6,61 7E-P2 

1, 103E-01 

7 , 29  St-02 

5, 457k -0  3 

6, 587 E -02 

tem 

7.2U0E-04 

4, 404E-04 

2nd 

2 

2.247F-05 

1.991E-05 

3,8856-05 

6,01 5E-05 

7,2686-05 

8,61  IE-05-1,. 1126-04 

5.260E-04 

l,926€-»3 

4.867E-03 

6, 931  f *k,2 

7,1 3 IE-02 

3,5066-02 

1.294E-03 

6, 4496-02 

tem 

1, 346E-04 

l,  180C-04 

6, 3406-05 

6,8076-04 

3rd 

3 

8,6t6E-05 

7, 353F.-0S 

5,785E-05 

5.026E-05 

3.694E-05 

3, 86 l E-04 

1 , 38  3E“03 

6.317E-03 

a,577F*0S 

1.657E-02 

2,6«2t*02 

1,69  76-03 

l , 900E-03 

tem 

3,6046-04- 

>2,7046-04 

6,5001-05 

3,  746E-04 

lith 

4 

2,6636-06 

4.364E-06 

l,57u€-07 

4,8756-85 

1,71 2E*04 

2,9126-04 

1,0*46-03 

2,9676-03 

l .519E-03 

2, 726E-0J 

2, 744E*03 

9 , 66  3E-04 

6 , 8031-04 

tem 

2,4436-04 

2,6306-04 

5,2fc56-*5 

1 , 1891-04 

5th 

5 

1.092F-06 

4 , 884f -06 

3, 95 Jt-06 

4 , 1 89E-05 

4,590E-05 

1,421 t -04 

2,95bE-04 

t ,4566-03 

5,2*26-04 

6 , 337E*04 

1,0*16-03 

3,72/6 *04 

7 , 95UF -05 

tem 

1,0491-04 

2.04RE-04 

9, 805  F.  *06 

1,30*6-05 

6th 

6 

3, 1 37F-06 

l ,928f *06 

2. J59F-06 

1.429E-05 

4,61 9E*0b 

2,8336*05 

5.238F-05 

1 ,6256-04 

l,93lf*04 

1 , 868E "04 

1,2*96-04 

1 , 1 7lE»0a 

2.024L-05 

tem 

3, iuaE-05 

I ,11696-04 

Figure  B2.  Transmittal  Card  Format  for  Joint  Fits  (Annotated) 
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AFWL-TR-78-4 


APPENDIX  C 

FORMAT  OF  PRINTED  OUTPUT  FROM  FPSPFT 


The  printout  from  the  fitting  program  FPSPFT  is  organized  into  five  ta- 
bles, with  3ome  information  displayed  in  more  chan  one  table.  Samples  of 
these  tables  for  a 9-group  beta  spectrum  (chosen  to  limit  the  fifth  table  to 
a single  page)  are  shown  in  figs.  Cl  through  C5.  Again,  these  figures  are 
liberally  annotated  and  should  be  nearly  self-explanatory.  All  except  the 
fourth  table  use  the  same  three  lines  for  a page  heading. 

The  original  purpose  of  the  first  output  table  was  to  summarize  the  strug- 

2 

gle  to  find  the  absolute  minimum  in  x » and  that  is  what  we  illustrated  in  fig. 
Cl.  In  practice,  however,  the  listings  delivered  to  AFWL  are  the  product  of 
the  additional  analysis  described  in  Sections  VI. 3 and  vi.6,  and  hence  are 
made  by  dictating  fixed  decay  constant:-  and  refitting  only  the  amplitudes. 

.cordlngiy,  they  include  further  information  between  lines  3 and  4 of  fig. 

Cl  to  identify  the  imposed  set  of  decay  constants,  and  the  table  contains 
only  the  first  of  the  six  columns. 

The  second  output  table,  as  shown  iu  fig.  C2,  displays  the  fitted  param- 
eters. For  the  free  fits,  the  transmittal  cards  are  essentially  equivalent 
to  reading  this  table  across  the  rows  in  turn.  For  the  joint  fits,  the  col- 
umns headed  lambda  are  punched  first  (they  are,  of  course,  the  same  for  each 
energy- isotope  case).  Then  the  amplitudes  are  punched,  again  reading  down 
the  columns. 

The  third  table  contains  the  most  detailed  information  regarding  the  fits, 
and  hence  is  too  crowded  to  explain  each  column  directly  in  fig.  C3»  The 
meanings  of  the  column  headings  are  fairly  evident,  but  require  some  additional 
explanation: 

GROUP  Number  of  the  group. 

MEV  Upper  energy  of  the  group. 

ITER  Total  number  of  iterations  to  reach  convergence.  In  the  final  list- 

ings this  appears  as  0 because  the  decay  constants  were  dictated. 

DAMP  Smallest  value  of  the  damping  constant  (see  equation  A7)  required  to 

achieve  the  fit  to  this  group.  Here,  again,  the  final  listings  al- 
ways give  1.0  since  no  iteration  was  required.  Below  the  value  of 
DAMP  for  each  group,  this  column  is  used  to  number  the  times  in  the 
SECONDS  column. 
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Figure  C3.  Example  of  Beginning  of  Third  Output  Table  (Annotated) . The  column  headings 
are  defined  in  this  appendix. 
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DATA  Intensities  in  this  group  calculated  by  CINDER/FPSPEC.  These  are 

in  particles/ (MeV1 3* fission) , not  in  particles/ (s* fission) . 

PCT.DIF.  FIT-DATA  expressed  in  percent  of  DATA.  The  rms  value  of  PCT.DIF. 

I . for  this  group  appears  below  the  last  time-step. 

i FIT  Intensities  calculated  from  the  fitted  parameters,  in  the  same  units 

as  DATA. 

TAU  Empirical  mean  life  (reciprocal  of  the  decay  constant)  in  seconds, 

determined  by  fitting  a Lagrange  polynomial  to  the  logarithm  of  the 
data,  as  described  in  Section  VI. 1.  For  sparse  data  sets  the  ap- 
parent mean  life  implied  by  pairs  of  adjacent  times  is  used  instead.* 

CURV  Normalized  first  derivative  of  TAU. 

DEP  Product  of  CURV  and  the  normalized  2nd  derivative  of  TAU. 

LAMBDA  Fitted  decay  constants,  in  units  of  s-^-. 

AMPLITUDE  Fitted  amplitudes,  in  particles/ (s* fission) . 

FRACT.  Fraction  of  the  integrated  yield  from  0 to  infinite  time  that  comes 
from  this  term. 

Between  the  columns  labelled  FIT  and  TAU  appear  asterisks  and  numerals. 
These  are  intended  to  be  read  as  graphical  marks  against  the  time  scale  implied 
by  TAU.  The  numerals  mark  the  assumed  mean  lives  with  which  the  fit  began, 
and  the  asterisks  show  how  far  these  migrated  during  iteration.  Both  are 
rounded  to  the  nearest  row  in  the  TAU  scale. 

' The  fourth  output  table  (fig.  C4)  shows  the  comparison  between  the  sum 

of  DATn  over  all  groups  and  the  sum  of  FIT  over  all  groups.  The  three  types 
of  rms  error  are  discussed  in  Section  VI. 4. 

The  final  table  is  illustrated  in  fig.  C5.  All  of  the  information  in  this 
table  is  reconstructed  from  the  fitted  parameters.  In  principal,  almost  all 
of  it  occurs  in  the  third  table  also,  but  here  we  have  converted  the  units  to 
group  constants  by  multiplying  each  fitted  intensity  by  the  width  of  its  group. 
We  have  also  extrapolated  the  table  backwards  to  zero  time  and  forwards  to  10^ 
seconds.  We  have  discussed  the  accuracy  of  this  extrapolation  in  Section  VI. 5. 
The  total  number  of  particles  per  fission  integrated  to  infinite  time  appears 
only  in  this  fifth  table. 

* 

A rule  of  thumb  that  was  widely  used  in  the  very  early  days  of  Civil  Defense 

states  that  the  current  decay  constant  is  equal  to  the  reciprocal  of  the  time 

: since  fission.  This  approximation  implies  a t"-*-  decay  curve.  The  difference 

I between  SECONDS  and  TAU  is  a measure  of  the  inaccuracy  of  the  t”l  approxima- 

| tion. 
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Figure  C4.  Example  of  Fourth  Output  Table  (Annotated).  This  table  compares 
the  sum  of  the  fitted  spectrum  to  the  sum  given  by  CINDER. 
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